Improved Measurement of the Hydrogen 1S — 2S Transition Frequency 
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We have measured the IS — 2S transition frequency in atomic hydrogen via two photon spec- 
troscopy on a 5.8 K atomic beam. We obtain /is-2S = 2 466 061 413 187035 (10) Hz for the hyperfine 
centroid. This is a fractional frequency uncertainty of 4.2 x 10 _lj improving the previous measure- 
ment by our own group [M. Fischer et a/., Phys. Rev. Lett. 92, 230802 (2004)] by a factor of 3.3. 
The probe laser frequency was phase coherently linked to the mobile cesium fountain clock FOM 
via a frequency comb. 

PACS numbers: 06.20.Jr, 32.30.Jc, 42.62.Fi 
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For the last five decades spectroscopy on atomic hydro- 
gen along with its calculable atomic structure have been 
fueling the development and testing of quantum electro- 
dynamics (QED) and has lead to a precise determination 
of the Rydberg constant and the proton charge radius |l| . 
The absolute frequency of the 1S—2S transition has been 
measured with particularly high precision, so that it now 
serves as a corner stone in the least squares adjustment of 
the fundamental constants [2J. The resonance has been 
used to set limits on a possible variation of fundamental 
constants [3| and violation of Lorentz boost invariance [4j . 
It further promises a stringent test of the charge conjuga- 
tion/parity/time reversal (CPT) theorem by comparison 
with the same transition in antihydrogen [5j. 

In this Letter we present a more than three times more 
accurate measurement of the IS — 2S transition as com- 
pared to the previous best measurements [3J, Iff reaching 
the 4 x 10 regime. The key improvements are the 
replacement of a dye laser by a diode laser system with 
improved frequency stability for the two photon spec- 
troscopy. In addition, a direct measurement of the 25* ve- 
locity distribution of the thermal atomic hydrogen beam 
allows a more accurate characterization of the second or- 
der Dopplcr effect. Finally, a quench laser resetting the 
population to the ground state right after the hydrogen 
nozzle was introduced. This removes possible frequency 
shifts due to the high density of atoms and a possible dc 
Stark shift within the nozzle. 

The extended-cavity diode spectroscopy laser system 
(ECDL) has been described elsewhere [7|. An ECDL 
master oscillator near 972 nm is amplified in a tapered 
amplifier. The laser radiation is frequency doubled twice 
within two resonant cavities to obtain 13 mW of the re- 
quired uv light near 243 nm. Locking the laser to a high 
finesse cavity made from ultra-low expansion glass leads 
to a line width of less than 1 Hz and a fractional fre- 
quency drift of 1.6 x 10~ 16 s _1 [8]. A fiber laser frequency 
comb with a repetition rate of 250 MHz is used to phase- 
coherently link the cavity frequency to an active hydro- 
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FIG. 1: (color online). Schematic of the beam apparatus. A 
standing laser wave at 243 nm (between grey mirrors) with 
a 1/e 2 waist radius of wo = 292 (tm at the flat cavity front 
mirror excites the sharp 15* — 2S transition in a co-linearly 
propagating cold thermal beam of atomic hydrogen emerging 
from a cooled copper nozzle. The 25 state is detected af- 
ter quenching with a localized electric field which releases a 
Lyman-a photon. See text for further details. 



gen maser which is calibrated using the mobile cesium 
fountain atomic clock FOM [10J. We apply cycle slip 
detection as described in j9| ■ 

The beam apparatus (Figjp follows the design of [11| 
and has been used before [12|. Hydrogen is dissociated 
at a pressure of 1 mbar in a radio frequency (rf ) dis- 
charge running in a sapphire tube. A Teflon capillary 
controls the flow and Teflon tubing guides the atomic 
hydrogen to a copper nozzle cooled to 5.8 K by a liquid 
helium flow cryostat. The dissociation fraction at the 
nozzle (2.2 mm diameter) is about 10%. After 45min 
of operation the nozzle closes up with frozen molecular 
hydrogen. We then heat it to 20 K for 15min to evapo- 
rate the hydrogen molecules. The atomic beam is defined 
by a 2.4 mm (front) and a 2.1mm (rear) diameter aper- 
ture which also separate the differentially pumped exci- 
tation region (10 _5 mbar / 10 _8 mbar). This region is 
shielded by a grounded Faraday cage from electric stray 
fields that may build up from laser ionized hydrogen. An 



enhancement cavity with a finesse of 120 forms a stand- 
ing wave for collinear excitation of the 15 — 25 transition 
by two counter-propagating photons. For detection, ex- 
cited, metastable atoms are quenched via the 2P state by 
an electric field (lOV/cm) and the emitted 121 nm pho- 
tons are detected by a photon multiplier tube (PMT). 
The intra-cavity power is monitored by measuring the 
transmission of the cavity using a photo diode connected 
to an integrating sphere. 

During twelve consecutive days starting on May 30 
2010 with a break on June 9, 1587 15 — 25 spectra have 
been recorded. For each 15* — 25 line the spectroscopy 
laser's frequency samples the transition in random order. 
At each frequency point we count Lyman-a photons for 
one second (with a 50 % duty cycle due to the chopper) 
at two different 243 nm laser intensities. A double pass 
AOM in zcroth order placed in front of the enhancement 
cavity allows to quickly alter the power level under oth- 
erwise identical conditions. After switching the AOM 
we wait several milliseconds to avoid possible frequency 
chirps. 

We apply an external magnetic field of 0.5 mT to sep- 
arate the hyperfine components. For the spectroscopy, 
we use the transitions from F = l,mp = ±1 to F 1 = 
l,mpi = ±1 whose Zeeman shifts cancel to zero. The 
hyperfine centroid frequency can be obtained by adding 
A/hfs = + 310 712 229.4 (1.7) Hz calculated from the ex- 
perimental results for the 15 and 25 hyperfine split- 
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tings [13, 

Now, we discuss the compensation of the two main 
systematic effects and estimations of the remainder. The 
Doppler effect, due to the velocity v of the atoms is can- 
celled to first order by virtue of the two photon excita- 
15|. The remaining second order Doppler 
~ w2 /is-2s/(2c 2 ) is compensated in two 



tion scheme 
shift A / d p 



steps: First, we chop the excitation light at 160 Hz (see 
Fig.QJ which allows time of flight resolved detection. Af- 
ter the light is switched off we only start recording the 
25 signal after certain delays r to let the fastest atoms 
escape. This samples the slow tail of the velocity dis- 
tribution. The 25 counts are then sorted into twelve 
time bins n = 10... 210 /is, t 2 = 210 . . AIO/js, ..., 
Ti2 = 2210 ... 2410 jUS. This removes most of the sec- 
ond order Doppler effect, but is not sufficient to meet 
the current level of statistical uncertainty. Therefore, we 
measure the velocity distribution of the 25 atoms at the 
detection point as described below and use it to further 
correct the second order Doppler effect. The residual un- 
certainty of this procedure was determined by evaluating 
a second data set that was generated using a Monte Carlo 
(MC) simulation [16], [l_7j in exactly the same way. For 



various simulation parameters such as temperature, ge- 
ometry and initial 15 velocity distributions we find the 
uncertainty to be smaller than 2.0 x 10~ 15 which is be- 
low the current statistical uncertainty of 2.6 x 10 -15 (see 
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FIG. 2: (color online), (a) Low velocity part of the Doppler 
profiles of 2S atoms recorded for different delays r at 370 mW 
intra-cavity power, (b) Fractional second order Doppler cor- 
rection Afd P {P)/fis-2S versus intra-cavity power for delay 
re = 1010 /is . . . 1210 /ts. Each point represents the Doppler 
correction as calculated from the central velocity of a single 
velocity profile measurement along with a linear fit. 
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The second systematic effect that needs to be corrected 
for is the ac Stark shift which is mostly linear in laser 
power P. However, a small quadratic contribution [l6| 
must be taken into account, before we can apply a linear 
extrapolation using the stable, but otherwise not pre- 
cisely calibrated laser power readings. The main contri- 
bution to this non-linear ac Stark shift is due to ionization 
of the 25 atoms by a third 243 nm photon that removes 
preferably atoms that see larger laser powers. For excita- 
tion laser powers of 300 mW as present in our experiment 
the quadratic ac Stark shift contributes on the order of 
1 x 10 -14 as derived from the MC simulations. This is suf- 
ficiently small to rely on these simulations that assume a 
Maxwell distribution for the 15 atoms using the absolute 
laser power within 20 % relative uncertainty. Modeling 
and subtracting the delay dependent quadratic ac Stark 
effect in this way then allows to linearly extrapolate the 
line centers, without knowing the exact laser power cali- 
bration. This procedure reduces the overall ac Stark shift 
uncertainty to 0.8 x 10" 15 . 

We find the experimental (and simulated) line centers 
by fitting Lorentzians which represents a good approx- 
imation of the line shape for delays T4 = 610 . . .810 /is 
and higher. For lower delays the second order Doppler 



effect causes an asymmetry so we do not evaluate them. 
A small residual asymmetry for the longer delays is de- 
termined and taken into account by comparing with the 
MC simulation. Again, we use this simulation only for 
small corrections. 

To correct the second order Doppler effect A/d p = 
—v 2 fis-2s/{2c 2 ), an accurate understanding of the ve- 
locity distribution is desirable. Previously, this informa- 
tion has been extracted from the line shape of the 25 
spectra with an uncertainty of 8 x 10~ 15 |6||. Here, we 
measure the 25 velocity distribution directly via the first 
order Doppler effect on the 25 — AP one photon transi- 
tion which we excite under an angle of 45° near 486 nm 
(Doppler laser in Fig. [1]). The 25 — AP transition has a 
sufficiently narrow natural line width of 13 MHz (corre- 
sponding to Av = 8m/s at 45°) to resolve velocities on 
the level of 1 m/s. 

The 4P state decays to the ground state with a 90 % 
branching ratio emitting a 97 nm Lyman-7 photon which 
can be easily detected using a channeltron. Pulsing the 
486 nm Doppler laser with an AOM avoids power broad- 
ening (and loss of velocity resolution) while providing 
equal quench probability for atoms of different veloc- 
ity. We cross the 486 nm beam with the atomic beam 
right between the quench electrodes which are grounded 
for the velocity distribution measurements. Using the 
same delayed detection as for the 15 — 2S spectra al- 
lows to extract the velocity distribution of 2S atoms 
that contribute to the signal with delay r. Working 
with the excited state directly gives the advantage of 
measuring the convolution of the velocity distribution 
with the excitation probability making the simulation of 
excitation dynamics unnecessary for this purpose. The 
low velocity part of a typical Doppler profile for delays 
T4 = 610 . . .810 /is to T7 = 1210 . . . 1410 /us is shown in 
Fig. E] (a). 

From 131 recorded Doppler profiles p T (v,P) we cal- 
culate the second order Doppler effect for each delay r 



according to A/d p (P) 



:{P)fis-2s/(2c 2 ) where the 



central velocity v c is determined by a Gaussian fit to the 
velocity profiles (see Fig. [2](a)). The power dependence 
arises from ionization losses of slow atoms. A linear fit 
to A/d p (P) reveals the second order Doppler effect cor- 
rection as shown in Fig.[2](b). 

The uncertainty in the second order Doppler correc- 
tion is caused by three main sources. First, the statisti- 
cal uncertainty obtained from linear regression analysis 
of A/d p (P) contributes 1.7 X 10~ 15 . Second, during the 
velocity measurements the 1S—2S spectroscopy laser was 
kept on the resonance only within ±160 Hz. MC simu- 
lation reveals an associated uncertainty of 0.8 x 10 -15 . 
Third, the 45° angle between the atomic beam and the 
laser beam used to measure the velocity distribution can 
only be adjusted within ±1°. This translates to an uncer- 
tainty in the second order Doppler effect of 0.8 x 10 -15 . 
Summing in quadrature leads to an overall uncertainty 
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FIG. 3: (color online). The plot shows fis~2S with all 
corrections applied for delays T4 = 610 . . . 810 /is to ti = 
1210 . . . 1410 /is (filled points) along with the mean as cal- 
culated from delays T5,6 (see text). The error bars represent 
the \a statistical uncertainty. The open circles show the same 
data with the second order Doppler effect not corrected, the 
open squares do neither include the second oder Doppler ef- 
fect nor the quadratic ac Stark shift corrections. 



of the second order Doppler correction of 2.0 x 10 15 . 

The fully corrected data is shown in Fig. [3] The tran- 
sition frequency /is- 2s is independent of the delay for 
r 4 = 610 . . . 810 lis, to t-i = 1210 . . . 1410 /js unlike before 
the correction. The insignificant remaining slope can be 
readily explained within the uncertainty in the correction 
of the second order Doppler effect. For the final analysis 
we only use delays T§fi. For delays T1...4 the second order 
Doppler effect cannot be sufficiently characterized. For 
delay T7 and higher, the quadratic ac Stark shift correc- 
tion cannot be extracted with competitive uncertainty 
due to the inaccuracy of the absolute power measure- 
ment. Also, the statistics is poor for these delays. 

As mentioned above, the pressure shift and the dc 
Stark effect have been greatly reduced as compared 
to " 



12j by the introduction of the quench laser. With 



an atom flux of 10 17 particles per second we find the re- 
maining pressure shift within the excitation region to be 
well below 0.2 x 10~ 15 [ll|. In [T3| the remaining dc Stark 
shift was measured to be 0.4 x 10~ 15 . 

Scattered light from the 486 nm quench beam can cause 
an ac Stark shift on the 15 — 25 transition. We experi- 
mentally restrict such an effect on the 0.4 x 10 -15 level by 
intentionally increasing the stray light. Also, atoms leav- 
ing the quench beam experience a low intensity region at 
the far wing of the Gaussian beam profile at which they 
are shifted but not quenched. We numerically simulate 
this effect by extending our MC simulation by a spatially 
dependent quench rate and an ac Stark shift due to the 
quench laser. We find no significant shift with an unccr- 



tainty of 0.8 x KT 15 . 

The Zeeman shift of the F = l,mp = ±1 to f = 
X,mF' = ±1 hyperfine transitions is ±360Hz/mT and 
consequently averages to zero for equal populations in 
both mp components. In a dedicated experiment, we 
have measured the Zeeman shift to be below 0.38 x 10~ 15 
by varying the applied magnetic field. 

Following [20j|, we correct the shift due to blackbody 
radiation at 300(30) K as -4.1(1.2) x 10~ 16 . Measur- 
ing the temperature fluctuations of the power switching 
AOM with the fed rf power we find an upper limit for 
an associated frequency chirp at 1.1 x 10~ 16 . From a 
measurement of the power radiated from the rf discharge 
into the excitation region we can estimate a possible ac 
Stark effect below 1.2 x 10 -17 . Theoretical considerations 
limit the shift due to higher order modes in the enhance- 
ment cavity tol.2xl0~ 17 and that by line pulling by the 
mp = hyperfine component to 1.6 x 10~ 18 . A residual 
recoil shift can be restricted by analysis of the observed 
Doppler broadening to be below 3.6 x 10~ 18 . 

The uncertainty of the fountain clock has been eval- 
uated to 0.81 x 10 -15 [2l| . A conservative estimate of 
the height difference between clock and experiment con- 
tributes a gravitational red shift uncertainty of 7.7 x 
10~ 17 . The uncertainty budget is summarized in Tab. HI 

Summarizing all corrections and uncertainties we find 
the IS — 2S hyperfine centroid frequency to be fis-2S = 



TABLE I: Uncertainty budget, a - uncertainty. 





a 


c/flS-2S 




[Hz] 


[10~ 15 ] 


statistics 


6.3 


2.6 


2nd order Doppler effect 


5.1 


2.0 


line shape model 


5.0 


2.0 


quadratic ac Stark shift (243 nm) 


2.0 


0.8 


ac Stark shift, 486 nm quench light 


2.0 


0.8 


hyperfine correction 


1.7 


0.69 


dc Stark effect 


1.0 


0.4 


ac Stark shift, 486 nm scattered light 


1.0 


0.4 


Zeeman shift 


0.93 


0.38 


pressure shift 


0.5 


0.2 


blackbody radiation shift 


0.3 


0.12 


power modulation AOM chirp 


0.3 


0.11 


rf discharge ac Stark shift 


0.03 


0.012 


higher order modes 


0.03 


0.012 


line pulling by uif = component 


0.004 


0.0016 


recoil shift 


0.009 


0.0036 


FOM 


2.0 


0.81 


gravitational red shift 


0.04 


0.077 


total 


10.4 


4.2 
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FIG. 4: (color online). The IS — 2S centroid frequency of 
the current measurement is compared with the previous two 
measurements fj, [g. The black points represent an average 
per day, the red points (with label) are the final value with 
the error bar representing the la uncertainty. 



2 466 061 413 187 035 (10) Hz. This corresponds to a frac- 
tional frequency uncertainty of 4.2 x 10 -15 and is in good 
agreement with our previous best measurement [3j but 
3.3 times more accurate. A comparison of the previous 
two measurements along with the current result is pre- 
sented in Fig.|4j The excessive day to day scatter present 
in the 1999 and 2003 measurements has been attributed 
to a dye laser instability and has consequently been re- 
moved in the current measurement. 
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